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Abstract 

We recalculate the leading one-loop contributions to fj, — ► ery and u — ► eee in the 
Littlest Higgs model with T-parity, recovering previous results for the former. When 
all the Goldstone interactions are taken into account, the latter is also ultraviolet 
finite. The present experimental limits on these processes require a somewhat heavy 
effective scale ~ 2.5 TeV, or the flavour alignment of the Yukawa couplings of light 
and heavy leptons at the ~ 10% level, or the splitting of heavy lepton masses to a 
similar precision. Present limits on r decays set no bounds on the corresponding 
parameters involving the r lepton. 



Contents 



1 Introduction 

2 The Littlest Higgs model with T-parity 

2.1 The Lagrangian 

2.2 Flavour mixing 

3 New contributions to LFV processes 

3.1 \i — > ej 

3.2 \i — > eee 

4 Numerical results 

5 Conclusions 
A Physical fields 

B Feynman rules 

B.l SM with massive neutrinos 

B. 2 LHT model 

C Loop integrals 

C. I Two-point functions 

C.2 Three-point functions 

C.3 Four-point functions 



1 Introduction 



Little Higgs models [1] offer an explanation to the little hierarchy between the Higgs mass 
Mh assumed to be near the electroweak scale v = 246 GeV and the new physics (NP) scale 
/, whose natural value is expected to be ~ 1 TeV [2]. In contrast with supersymmetry, 
where the large one-loop Standard Model (SM) contributions to the Higgs mass are can- 
celled by the contributions from the corresponding supersymmetric partners with masses 
~ 1 TeV and spins differing by ±1/2 (see [3] and references therein), Little Higgs (LH) 
models stabilize by making the Higgs a pseudo-Goldstone boson of a broken global 
symmetry. The cancellation is in this case between particles with the same spin belonging 
to the same multiplets of this approximate symmetry. Which of these SM extensions, if 
any, is at work will be hopefully established at the LHC [4,5]. 

Supersymmetry is linearly realized in the minimal supersymmetric SM (MSSM). This 
and other simple supersymmetric extensions of the SM have other interesting phenomeno- 
logical properties, like, for instance, the unification of gauge couplings at very high scales [6] , 
which is not the case for LH models. However, they have no built-in low energy mechanism 
to explain the observed fermion mass hierarchy or flavour conservation. As a matter of 
fact, it can be argued that supergravity is phenomenologically relevant [3,7] because it can 
provide the necessary initial conditions to explain the precise fine-tuning required among 
the many new parameters of the MSSM, which otherwise would result in too large flavour 
changing processes [8] . This has been historically the problem of many SM extensions [9] . 
If the NP is near the TeV scale, it faces in general the problem of naturally explaining why 
it is aligned with the SM Yukawa interactions, as experimentally required. Although the 
SM can not explain the large hierarchy between fermion masses, which by the way is several 
orders of magnitude more demanding than the little hierarchy, it naturally accommodates 
the absence of flavour changing neutral currents (FCNC) [10]. LH models are not designed 
to solve the flavour puzzle either, and one must expect stringent constraints on the new 
parameters involving the heavy sector. The study of FCNC processes in Littlest Higgs 
models has been addressed in the literature [11]. In this paper we revise the calculation 
of the decay rates of the lepton flavour violating (LFV) processes /x — > [12, 13] and 
/i — > eee [13] in the Littlest Higgs model with T-parity (LHT) [14], obtaining an ultravio- 
let finite result also for the latterQ Indeed, when all Goldstone boson interactions of the 
new leptons are taken into account, the one-loop contributions to the amplitudes are well- 
defined [16,17], scaling approximately in the two family case like (t> 2 // 2 ) sin2# 5, where 
6 is a measure of the misalignment between the heavy and SM lepton Yukawa couplings 
and 5 is the corresponding heavy lepton mass splitting. As a consequence, the present 
experimental limits require fine tuning the Yukawa couplings of the new heavy leptons up 
to 10%, aligning them with their SM counterparts, or making the heavy masses quasi- 
degenerate. One might also rise the NP scale degrading the motivation of the LH scenario 

^ee also [15] which appeared when we were preparing this manuscript. There the cancellation of 
ultraviolet divergences in the LHT model is also shown, flavour violation in the quark sector is explored 
and the phenomenology of K — > kvv is analyzed. 
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itself. In the general case with three families the new contributions must be tuned to a 
similar precision but the parameter dependence is more involved. The calculation also 
applies to r decays, but the corresponding limits are not restrictive at present. Moreover, 
it can be easily extended to \i — e conversion in nuclei [18]. A complete phenomenological 
analysis comparing as well different LH models will be presented elsewhere. 

In LH models the Higgs is a pseudo-Goldstone boson. Thus, Mh is naturally small as 
long as the new scale / is relatively low, because one expects that cancellations are only pro- 
tected to one loop and for the dominant contributions. Hence, Anf can not be much larger 
than 10 TeV if we do not wish to invoke some fine tuning again. However, as the model 
introduces heavy particles the new one-loop contributions to electroweak observables may 
require rising / significantly above 1 TeV in the absence of model dependent cancellations, 
in order to be consistent with present electroweak precision data (EWPD) [19]. The LHT 
is an economical realization of the LH scenario with the further virtue of keeping the new 
contributions to EWPD small. It incorporates a discrete symmetry under which the new 
particles are odd and the SM ones even. Then all vertices must have an even number of new 
particles, if any. Similarly to the R symmetry in supersymmetric models, the T symmetry 
allows us to weaken the experimental limit on the LH effective scale below the TeV [20]. 
This symmetry also makes stable the lightest T-odd particle, offering, like R-parity does 
in the supersymmetric case, an alternative candidate for cold dark matter [21]. 

Nevertheless, as already emphasized these models are not a priori designed to deal 
with the flavour problem. Therefore, it is important to investigate the constraints on the 
model parameters implied by the stringent experimental limits on FCNC. We follow an 
operational approach and calculate the leading contributions to /i — > and \i — > eee 
in the LHT, recovering previous results for the former [12, 13] but an ultraviolet finite 
result for the latter. We focus on these processes because the lepton sector is free from 
large strong corrections, and the experimental limits are quite demanding. In Section [2] 
we review the LHT model to introduce the notation and the Feynman rules needed. The 
one-loop amplitudes of the LFV processes \i — > e'j and // — > eee are discussed in Section [31 
The calculation is straightforward but cumbersome, requiring a careful bookkeeping of the 
different terms. In Section @] we present the numerical results discussing the dependence on 
the different parameters of the model. Finally, Section [5] is devoted to conclusions, where 
we also briefly comment on r decays. 

2 The Littlest Higgs model with T-parity 

2.1 The Lagrangian 

The LHT is a non-linear a model based on the coset space SU(5)/SO(5), with the SU(5) 
global symmetry broken by the vacuum expectation value (VEV) of a 5 x 5 symmetric 
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The 10 unbroken generators T a , which leave invariant S and then satisfy T a E +S (T a ) T 
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diag(2,2,2,-3,-3), 
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with a a the three Pauli matrices. The VEV in Eq. (12.11) breaks this gauge group down to 
the SM gauge group SU(2) i xU(l)y, generated by the combinations {Q\ + Q 2 , Y\ + Y 2 } C 
{T a }. The orthogonal combinations are a subset of the broken generators, {Q± — Q 2 , Y\ — 
Y 2 } C {X a }. Thus, the Goldstone fields 
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decompose into the SM Higgs doublet (— i7r + /v^2, (v + h + i7r°)/2) T , a complex SU(2)x 
triplet and the longitudinal modes of the heavy gauge fields ^,00° and rj 
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2 In the following we use for the SM fields and couplings the conventions in Ref. [22]. In particular, 
b+ = -m+, 0° =7T°. 
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As emphasized in the previous section, we can make the new contributions to elec- 
troweak precision observables small enough introducing a T-parity under which the SM 
particles are even and the new particles are odd. An obvious choice for the action of such 
T-parity on the gauge fields Gi is the exchange of the gauge subgroups [SU(2)xU(l)]i and 
[SU(2)xU(l)] 2 , 



Go. 



(2.6) 



Then, T invariance requires that the gauge couplings associated to both factors are equal. 
This leaves the following gauge Lagrangian unchanged, 
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(2.7) 



where 



Bj^u = d^Bju — d u Bj^. (2.8) 



(Summation over index a, which runs on the corresponding SU(2) generators, is always 
assumed when repeated.) The T-even combinations multiplying the unbroken gauge gen- 
erators correspond to the SM gauge bosons, 
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expand the heavy gauge sector. 

In order to ensure that the SM Higgs doublet is T-even and the remaining Goldstone 
fields are T-odd, the T action on the scalar fields is defined as follows, 

n -nun, n = diag(-i,-i,i,-i,-i), (2.11) 

where f2 is an element of the center of the gauge groupjl which commutes with E but not 
with the full global symmetry. Then, 



and the scalar Lagrangian 
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Note that we have reversed the sign of f2 as compared to the literature, to make it a group element. 



6 



with 



(2.14) 



is also gauge and T-invariant. 

This discrete symmetry must be implemented in the fermion sector too. This is less 
straightforward. In fact, there is no proposed model fulfilling the three desired conditions: 
to give masses to all (SM) fermions with Yukawa couplings, preserving a discrete sym- 
metry under which all new particles are odd and the SM ones even, and keeping the full 
global symmetry before introducing the symmetry breaking. Although terms explicitly 
breaking the global symmetries at the Lagrangian level must manifest as badly behaved 
contributions to physical processes [17], this will not be our case since all the explicit 
couplings entering in the calculation we are interested in can be derived from Lagrangian 
terms which are symmetric. Following Refs. [23,24] we introduce two left-handed fermion 
doublets in incomplete SU(5) multiplets, one transforming just under SU(2)i and the other 
under SU(2) 2 , for each SM left-handed lepton doublet: 
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under an SU(5) transformation V. We define the T-parity action on these fermions 



(2.17) 



Then the usual T-even combination ^ + r2S \I / 2 remains light and is identified, up to the 
proper normalization, with the SM fermion doublet. The T-odd combination ^ — QH ^/ 2 
pairs with a right-handed doublet (eigenvector of T), in a complete SO(5) multiplet, 
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where U is an SO(5) transformation defined below, to form a heavy Dirac doublet. With 
this aim in mind, a non-linear Yukawa Lagrangian is introduced, 



C Yh = -«/ + *iS e f ) ^r + h.c. , 
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where £ = e lU ^ . This is indeed T-invariant, since Eq. (12.111) implies 



and invariant under global transformations, 
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where V is the global SU(5) tranformation and U a function of V and II taking values 
in the Lie algebra of the unbroken SO(5). It must be noted that the gauge singlet xr, 
completing the SO (5) representation 
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and assumed to be heavy, is T-evenQ On the other hand, the extra doublet xJjr, which i 
also assumed to be heavy enough to agree with EWPD, is T-odd as desired. 

We have just introduced all heavy fields we need. However, one important comment 
is in order. The Yukawa-type Lagrangian Cy H fixes the transformation properties of the 
heavy fermions and then their gauge couplings, in particular the non-linear couplings of 
the right-handed heavy fermions [25], 
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The Lagrangian of Eq. f!2.23|) includes the proper 0{v 2 / f 2 ) couplings to Goldstone fields, 
absent in [13,26], that render the one-loop amplitudes ultraviolet finite. Besides, in order 
to assign the proper SM hypercharge y — — 1 to the charged right-handed leptons £r, which 
are SU(5) singlets and T-even, one can enlarge SU(5) with two extra U(l) groups, since 
otherwise their hypercharge would be zero. Then, the corresponding gauge and T invariant 
Lagrangian reads 



C' F = ii R r(d^ + ig>yBje R . 



[2.25) 



4 If we had defined the T action on the fermions *S?i < — ► — ^0^2, ^r. — ► —^r and the Yukawa 
Lagrangian with Si's, Cy H = —re/ (^C + , 3 > iSo^^^) &r + h.c., all new fermions would be T-odd and the 
new Lagrangian invariant under the new T-parity [24] , but not under the full global symmetry because ft 
does not commute with SU(5) neither with SO (5), although it does commute with the gauge group. We 
must insist that the explicit couplings entering in our calculation are the same in both cases. 
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For the lepton sector and the calculation we are interested in these are all the necessary 
Lagrangian terms. However, in order to define what a muon or an electron is, we have to 
diagonalise the mass matrix (Mt)ij = {^e)ijV in the corresponding Yukawa Lagrangian Cy 
which we assume to have all required properties [25,27] (and also to include light neutrino 

masses). This gives to leading order a mass term for the charged leptons rri(iT L t R + h.c, 
with 

mM' = CV?)i'«(A*)«(v£)tf « (2-26) 

and Vi R two unitary matrices^ 

Finally, in order to perform the calculation in the mass eigenstate basis we have to 
diagonalise the full Lagrangian 

C = C G + C s + C Yh + C F + C' F + C Y , (2.27) 

and reexpress it in the mass eigenstate basis. The corresponding masses and eigenvectors 
up to order v 2 /f 2 are given in Appendix^ The Feynman rules are collected in Appendix [Bl 
They are obtained expanding C to the required order. The coupling overlooked in [13] is 
the v 2 /f 2 correction to the right-handed coupling g R of the Zv l H v 3 H vertex, resulting from 
the expansion of the last two terms of Cf in Eq. (12.231) . 

2.2 Flavour mixing 

The new contributions to charged LFV processes must be proportional to the ratio of 
the electroweak and the LHT breaking scales v 2 /f 2 and to a combination of the matrix 
elements describing the misalignment of the heavy and charged lepton Yukawa couplings. 
Let us then set our conventions for the description of the heavy-light mixing relevant to 
our analysis, and in particular to the Feynman rules discussed above and collected in 

5 Right-handed leptons, as the other right-handed SM fermions, are usually taken to be singlets under 
the non-abelian symmetries, transforming only under the gauge abelian subgroup. We must note that 
this may be a too strong assumption. If we want to couple them to their left-handed counterpart, one 
may be inspired by the following observation. There is only one SU(5) singlet in the decomposition of 
the product of two E's and one left-handed fermion multiplet, YJ a =1 e a i a ^ a 3 a i a a [(S) QlQ2 (Yj) a . iai ^2a & + 
(S^) Q , lQ , 2 (S^) Q , 3Q , 4 4'i Q , 5 ], where £ a i a 2«3040 5 j g the totally antisymmetric tensor and the second term is 
the T transformed of the first one. Alternatively, one could multiply three E's and the other left-handed 
fermion multiplet, E„ 1= i e QlQ2Q3Q4Q5 ^ 6Q 1(S) aiQ2 (E) Q3Q4 (E) Q5a6 * lQ7 + (Et) QlQ2 (Et) Q3Q4 (Et) Q5Q6 * 2Q7 ], 
with 5 aea7 the Kronecker delta. In both cases, we get the wrong Higgs coupling. This is so 
because this product is an SU(5) singlet and then the Higgs coupling reads + (v + h + 
m°)y/\/2. (In these expressions there are neither f2's nor Eo's because the determinant of fl is 1 and 
e aiQ2Q3Q4Q5 (So)a 1 /3 1 (Eo) a2/ 3 2 (Eo)a3/33(So)a 4 /3 4 (Eo) Q5 /3 5 = fP^^K) Then, getting the correct coupling 
(v + h-iir°)l-/V2 + in v requires the explicit breaking of SU(5). If £ is introduced in the game, one 
eventually has to break SO(5) as well. 

6 We denote the mass eigenstates with primes when necessary to distinguish them from the current 
eigenstates. 
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Appendix [HI The SM interaction and mass eigenstates are related by the unitary matrices 
in Eq. (T2T26j) . 



e L = vte L , t R = v^ R - (2.28) 

Then the SM charged current Lagrangian reads 

Tc = ~^f^h + h.c. = -jf' L V^VlW^' L + h-c, (2.29) 

where we have also introduced the corresponding rotation for the neutrinos. Thus, only 
the combination Vp MNS = V^}V[ is observable. It must be noted, however, that the 
neutrino contributions to LFV processes are negligible in the SM because so are their 
masses. Hence, V£ can be assumed to be unity. Similarly we can also diagonalise the 
heavy Yukawa couplings in Eq. (12.191) . 

m i^i>i> = {VL%i K a( V R )jf Af> (2-30) 

where acts on the left-handed fields and acts on the right-handed fields. Note that 
there is no distinction between up- and down-type leptons. The T-odd gauge boson inter- 
actions arising from the corresponding kinetic terms for left-handed leptons in Eq. (12. 23f) 
are proportional to 

lL-(fU.L+ + h-c = ThlVFQh \ ^ | + li-c-- (2.31 ! 

where G- and II- are the heavy, T-odd gauge bosons and fermions and Il+ are the SM, 
T-even fermions in the interaction basis, whereas Gh = A H , Z h ,Wh] Ih = (v H ,£ H ) T ; 
and vl and are the corresponding mass eigenstates. Then, in analogy with the PMNS 
matrix, the observable rotations are now 



V Hv = V^Vl V m ee Vf X Vl (2.32) 

Note that both matrices are related, Vjj v Vhi = ^mns [28], but this relation can not 
be tested unless Vh v can be measured. The new contributions to the LFV amplitudes 
describing a muon decay to an electron are then proportional to V^V^, with i counting 
the heavy lepton doublets. 



3 New contributions to LFV processes 

As noted above, the SM contributions to the LFV processes fi — > ej and fi — > eee are 
negligible for they are proportional to the observed neutrino masses. On the other hand the 
new LHT contributions can be a priori large. In particular, one expects that the dominant 
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Figure 1: Effective vector-fermion vertex. 

contributions come from the exchange of the new vector bosons and heavy fermions required 
to realise the discrete symmetry T[] Here we study both processes in turn. 

The amplitude /i — > e7 is proportional to the vertex in Fig. [TJ whose most general 
structure for on-shell fermions fij can be written in terms of six form factors: 

il^Cpi.pa) = ie[^{FlP L + FlP R ) + {iFl + Fl l5 )a^Q v + {iF^ + Fj 75 )Q"], (3.1) 

with Prl — |(1 ± 75) and Q = p 2 — Pi the vector boson momentum entering into the 
vertex. If the vector boson V is a photon, the U(l) gauge symmetry is unbroken and 
current conservation implies 

{m i -m j ){Fl + F] i )+2iQ 2 F2 = 0, (3.2) 
-(m i + m j ){Fl-Fl)+2Q 2 F] 1 = 0. (3.3) 

Hence, the LFV process fj — > /j 7 with i 7^ j where the photon is on-shell (Q 2 = 0) 
is completely described by a dipole transition. Indeed, according to Eqs. (I3.2f3.3p F~l = 
Ffi = for on-shell photons, while the form factors Fg p do not contribute to the amplitude 
because real photons are transverse. Then, the total width for £j — > ^7 is given by 
[29, 32-35] 

r(^-^7) = fm|(|F2 / | 2 + |^| 2 ). (3.4) 

On the other hand, two types of diagrams contribute to fi — > eee (see Fig. [2]). Now, 
in the diagrams of the first type (penguins) the exchanged gauge boson V can be a 7 
or a Z but not a heavy vector boson for the coupling is forbidden by T-parity. (Higgs- 
penguins are neglected.) FY and FY do not vanish in penguins. In fact, for 7 these 
form factors are proportional to Q 2 , as we have explicitly checked. Besides, as the gauge 
boson couples to two on-shell electrons, the contributions from F^ P are irrelevant for they 
multiply the electron mass. The amplitude for this process also receives contributions from 
box diagrams. The total amplitude for fi{p) — > e(pi) e(p 2 ) ^(^3) can be then written [32] 

-'^7— penguin + M Z —penguin + A^box, (3.5) 

7 The addition of new vector- like leptons in general imply large FCNC already at tree level [10,29], and 
stringent constraints from EWPD and LFV processes [30]. In the LHT they are absent because T-parity 
forbids the coupling of a SM gauge boson to one light and one heavy fermion. Analogously, the presence 
of heavy scalar triplets with hypercharge 1 in general allows for their direct coupling to two (SM) lepton 
doublets (for a review and further references see [5,31]). This is also absent in the LHT because the triplet 
i> in Eq. (j2~5)) is T-odd and the SM leptons are T-even. 
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Figure 2: Generic penguin and box diagrams for \i — > eie 2 e 3 . Crossed diagrams with ei 
and e 2 exchanged must be added. 

with 

AVpcnguin = ^"(Pi) [Q 2 r(ALp L + A?P R )+m^Q v (A%P L + AfP R )]u(p) 

~xu{p 2 )l^v{p 3 ) - (Pi <-> P2), (3-6) 

-Mz-penguin = ^(Pl) {^{FlPl + F R P R )] U { P ) u(p 2 ) [i^Pl + Z R P R )\ v{p 3 ) 

z 

-(pi^Pa), (3-7) 

M hox = e 2 B{ [u{ Pl )^P L u{p)) [u{p 2 ) ltl P L v{p 3 )) 
+e 2 B* [u( Pi )YPru(p)} [u(p 2 )^P R v(p 3 )\ 
+e 2 P 2 L {{u{ Pi )YPlu{p)} [u(p 2 )j„P R v(p 3 )} - (pi <-> p 2 )} 
+e 2 P 2 R {[^(p^P^p)] [«(p 2 )7^(p 3 )] ~ (Pi <- P2)} 
+e 2 P 3 L {[iZ(pi)P L «(p)] [w(p 2 )P L t;(p3)] - (Pi <-> P 2 )} 
+e 2 P 3 R {[w(p!)P^(p)] [«(p 2 )P^(p 3 )] - (pi <-> p 2 )} 
+e 2 P 4 L {[u(p 1 )a^P L ?i(p)] bi(p 2 )a^P L t;(p 3 )] - (pi ^p 2 )} 
+e 2 Pf {[u(pi)^P^(p)] [«(p 2 )^P^(p 3 )] - (Pi <- P 2 )} • (3.8) 

We have defined new vertex form factors in the penguin amplitudes 

A{ = FHQ\ A? = FHQ\ A L 2 = (i£ + iP^)/m„ A* = (i£ - iP^)/m M , 

F L = -Ff, F R = -Ff, (3.9) 

and used that Q 2 <C M| in Eq. (13. 7p . zJI^ are the corresponding Z couplings to the 
electron in the SM (see Appendix [B]). The dipole form factors F^ E are dropped from the 
amplitude because their contributions are effectively suppressed by a factor m 2 /M^. The 
total width can then be written as [32,35]: 
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+ Sf* + A? B** + A L X B^ + Af B** + h.c.) 
- -(A$B%* + A^f* + A^f* + AfB%* + h.c.) 
+ ^{2(|F Li | 2 + \F RR \ 2 ) + \F LR \ 2 + \F RL \ 2 

+ {B{F* LL + B*F* RR + B L 2 Fl R + B*F* RL + h.c.) + 2(A\F* LL + Af F* RR + h.c. 
+ (A^Fl R + Af + h.c.) - 4(4?*^ + + h.c.) 

-2(A^F* RL + A^F* LR + h.c.)} 



(3.10) 



where 
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Note that the amplitude for the Z-penguin could have been cast into the box structure 
replacing 
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B« 



B{ + 2F LL 
Bf + 2Fj 



RR, 



B 2 + F LR , 
Bf + Frl- 



(3.12) 
(3.13) 
(3.14) 
(3.15) 



The branching ratios for both types of processes are obtained dividing by the SM decay 
width 
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G, 



na 



V2s 2 w M^ 



(3.16) 



For t decays the SM branching ratio must be corrected multiplying by 0.17 to take into 
account other possible decay channels. 



3.1 fi — > e7 

Let us now summarize the calculation for /i — > erf. The new one-loop Feynman diagrams 
contributing to the V/j,e vertex in the LHT model in the 't Hooft-Feynman gauge are listed 
in Fig. [3j They are classified in six topology classes. As explained above, in the decay 
fi — > e7 the photon is on-shell and then only the dipole form factors F]^ E contribute. 
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Figure 3: New one-loop diagrams contributing to V/xe in the LHT model. 
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They are proportional to the muon mass, reflecting the chirality flip character of the dipole 
transition. (The electron mass is neglected.) We separate the contributions exchanging 
Wh, Zh and Ah, expressing the results in terms of standard loop integrals (Appendix O). 



Diagrams exchanging Wh 



Taking Mi = M Wh and M 2 = vn v i^ and introducing the mass ratio 



m 



Hi 



(3.17) 



with niHi = Trip, 
W H (see Fig. Ejf 



m, 



we find the following contributions from diagrams exchanging 



II: F m\w h 



IV: 

V: 
VI: 

Total: 



T?l I 

r mWh 



-' lF E\W H 



-' iF e\w h 



aw 

Win 

16tt ^ 



a w 

Win 

16tt ^ 



£v&Vg [3CH-CJ 



Ev ie *v 



V y. 



m\w h 

7 



-iF7 



t i 



aw 
167r' 



77>7 I 



. 7 «h/ ra p 



(3.18) 

(3.19) 

(3.20) 
(3.21) 

(3.22) 



where aw = ot/s ^ and 



5 3a: — 15x 2 — 6x 3 



2d - 3Cn -x(C + 3C 1 + -C u 



6 



12(1 -x) 3 



+ 



3x 3 



2(1 -x) 



lnx. 



(3.23) 



The constant term drops from the amplitude due to the unitarity of the mixing matrix. 
This result is in agreement with [13,36,37]. 

It may be worth to note that these contributions are completely analogous to those of 
the SM with massive neutrinos, replacing Wh by W, u l H by z/j and Vhi by Vp MNS . For tiny 
neutrino masses, Xi = mljM w <C 1, 

5 



x 



F w (x)^--^ + 0(x 2 ), (3.24) 
and we recover a well known result [36] bounded by neutrino oscillation experiments: 



e7) 



3a 
32^ 



5>. 



PMNS^PMNS X i 



< 10 



-51 



(3.25) 
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Diagrams exchanging Zh 

Taking now Mi = Mz H and M 2 = rag, , with the same y^ we get: 

aw 



where 



F? 



m\z h 



F 7 



I: 

r M\Z H 

Total: Fli\z H 
F*(x) 



-iF? 



E\Z H 



-iF? 



E\Z H 



167T ^ 



32^ 



'Hi 



7 a w m M 



i 

i 

T, V mK e F z (y, 



Cq + 3Ci + 



He Vl 



M{ 



Co + 3Ci + — f^* 1 — 2 < ^ 11 



1 2x + 5x 2 — x 3 
3 + ' 



3x 2 



(1-x) 3 4(1 -x)' 



hire, 



in agreement with [13,37]. 



Diagrams exchanging Ah 

This contribution can be obtained from that of the diagrams with a Zjj, replacing 
Ajj. It is convenient to introduce the mass ratio 



Vi = ay%, a 



Then, 



T7~< I 
r M\A H 



M Wh 5c 2 



w 



M 2 s 2 
1V1 A H S W 



a w m„ 1 s 



16ttM 2 h 25 c 2 ,, 



a w m M l 
16V M^5 



Branching ratio 

Using M w /M Wh = v 2 /(4f 2 ) and M Wh = M Zh , we finally obtain: 



^ E ^ ( F ^) + F ^ + 



with Fw and Fz given in Eqs. (13.231) and (13.291) . respectively. 
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Bl B2 B3 B4 



Figure 4: Box diagrams for fj, — ► eie2e3 in the LHT model. Crossed diagrams with ei and 
e2 exchanged must be added. 

3.2 ji — > eee 

The self-energy diagrams do contribute to in this process, and must be included to 
calculate the penguin diagrams in /x — > eee . On the other hand, apart from the box 
diagrams, only 7- and Z-penguin diagrams contribute to /i — > eee in the LHT model. This 
is so because A# and Z# do not couple to two ordinary fermions, as required by T-parity 
conservation. (Higgs-penguins vanish in the limit of massless electrons, as do F^ in this 
limit too.) For the sake of brevity we present our results grouping together the Wh, Z h 
and A H contributions, but we distinguish among the 7- and Z-penguins in Fig. [3] and the 
boxes in Fig. HI 

The 7-penguin 

The form factors F]^ and F% have the same expressions (13. 22113. 28113. 31|) as for an on-shell 
photon, since terms of order Q 2 can be neglected. The contributions to F^, which are 
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proportional to Q 2 ~ m 2 as expected, are detailed below. 

For Mi = Mw H , M 2 = m v i and yi = m 2 Hi /M'^ H as before, the diagrams with Wh yield 



F l\w H - ^ Vm v m G w(Vi 



a w Q 2 



zZ V m v m G w(Vi)i (3-33) 



with 



Gw(x) = -\ + Si + 6C m + x Qfl; + C„„ - Mfc) - (2d + ±C U \ Q 2 

5 x(12 + x-7x 2 ) x 2 (12- 10x + x 2 ) , . . 

G - W = -18 + 24(1 -,)■ + 12(1-,)' ln " (M5) 

Relations (fCl20l and fj022l have been used in lECHjl . Note that owing to the unitarity of 
the mixing matrix the ^-independent terms in GV(x) drop out (including the ultraviolet 
divergence). The SM prediction is obtained by replacing Wh by W, v H by i/j and Vhi by 

v pmns- 

For M\ = Mz H {= M\v H ), M2 = m e i H ~ m v i and the same yi, the contribution of 
diagrams with Zh is 

F ^n = ^J2 V ^<G z ( yi ) (3.36) 



2 



% <5 



E^/V^G-W^), (3-37) 



with 



Gz(x) = (l + |) ^-i + lB 1 + Coo-|M?C ) - Qc + C 1 + i(2 + x)C 11 )g 2 



-G^(x) + C»(^), (3.38) 



s 1 x(18- 11a; -x 2 ) 4 - 16a; + 9x 2 , . . 

Gl' s = — + — ; ^ — - ; -;— lna;. 3.39 

z v ; 36 48(1 -x) 3 24(1 -x) 4 V ; 



The relation fl(J.22[) has been used in Eq. (13.3 

Finally, the contribution of diagrams with Ah is obtained from that of diagrams with 
Z H replacing Z H by A H , and y< by y[ = 5c 2 w yi/ s 2 w : 

F2\a h = E V2?Vg (3.40) 
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The Z-penguin 

The Z dipole form factors Fj^ E (which are chirality flipping and hence proportional to the 
muon mass) can be neglected as compared to F^ . This is in contrast with the 7-penguin, 

m 2 l /M WH ~ F2, to be compared with F[ ~ 1. 



< 



for which QFl IE (~ QF^ E ) ~ Q 2 /M$ 
This justifies to neglect Fj^ E in the Z-penguin (13.71) . 



Taking M x 
obtain: 



Mw H , M 2 = rriHi and yi = m 2 Hi / Myy , and using the unitarity of Vm we 



l \ w E 



aw 



8tt swCw 



— E V m V m{ ~ (~ + Bx + 6C 00 - y t M Wii C ) 



y iC 2 w (B 1 + 2C 00 ) 
V 



+ 2 1 + 



4 2 1 



Ci 



00 



2 

+ ^ [1 + 4(C 00 - Coo + M Wh (Co - 2C ))] 



aw 



— V V ie *V ifl \ -2c 2 f A - In Mw A 



8tt s w cw 



37T Sv^CvK 



(3.41) 



with 



6-x 2 + 3x 
= h — — hi 2, 



1 — x 



l-x) 



(3.42) 



and 



f l\a h 



aw 



8n swcw 



1 + 



T 2 



4 



+ 



— 1 
%' 



2svy 
5c w 



87T SiyC^y 



2c 



Wj 



4 2 1 



1 s 2 

25 Cyy 



2 



T 5 



5../ 



C, 



00 



C, 



00 



2 



v[ + 10 



^i ff Co 



0. 



(3.43) 



(3.44) 



Here xh is a constant defining the mixing between the heavy neutral gauge bosons and 
function of the gauge couplings (see Eq. ( 1A.4I) ). We observe that the only contribution to 
the Z-penguins comes from the diagrams with Wh, and it is proportional to v 2 /f 2 . The 
potentially dangerous ultraviolet divergences proportional to yi have cancelled thanks to 
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the proper v 2 / f 2 corrections to the ui^W^Z and Zv l HR y l HR couplings. The corrections to 
the latter were not included in [13]. 

For completeness, we give the prediction for the Z-penguin in the SM with light massive 
neutrinos. Although in the LHT the heavy leptons are vector-like and the Z boson couples 
to both chiralities, the final form of the vertex is the same. This is more easily seen in the 
unitary gauge, where the heavy modes contribution is only given by diagram I in Fig. [3] 
and is proportional to the v 2 / f 2 correction to the Zv % HR v % HR coupling (see [26] for further 
discussion). Taking M\ = M w , M 2 = m Vi and %i = m 2 ./M w , and using the unitarity of 
^pmns we obtain: 

f l\w = -^r y~] ^pmns^pmnsI ~ 2c w ( _ o + Bl + ® C °° ~ x i M w c o) 

+ |(l-2c 2 w ) (B l + 2C m ) 
- ~ + Bt + 2C 00 - ^M W (4C + Xl C ) 

a W 1 \ r T rei rrUi* \ a Jl ( A 1 



167T SwCw 



aw 



£ ^pmns<ns{ " (A, - In ^ - 1 

+ 2{B l + 2C 00 ) - XiM w (AC + XiC ) 

/J^PMNS^PMNS x iHw{ x i)i (3.45) 



167T SwCw 

which is, of course, finite and in agreement with Ref. [33] for Q 2 = 0. 



Box diagrams 

There are eight different classes of box diagrams grouped in types A and B in the LHT 
model (Fig. H]). In the limit of zero external momenta (all internal masses are much larger 
than the muon mass) all of them have the same form, being proportional to a scalar 
integral over the internal momentum q. Indeed, omitting the corresponding denominator 
(q 2 - m 2 m f(q 2 - M 2 J 2 , with G = W, Z or A, 

Al : (pi\ -fP L (-^ + m Hi )YP L \p) (p2\l,P L {-<i + m Hi ) lfl P L \p 3 ) 

= j(pi\ rrrp L b> (p 2 | i.i^pl Ip 3 > 

= 9 2 (pi|7^b)<P2|7^b3>, (3.46) 
A2 : - (pi| ^P L {-4 + m m )P L \p) (p 2 | P R {~4 + m Hi )^P L \p 3 ) 

= -m Hi m Hj (pi| 7 M P L \p) (p 2 \ I^Pl |ps) , (3-47) 
A3 : - {pi\P R {-4 + m m )^P L \p) (p 2 \ 1»Pl{-4 + m Hi )P L \p 3 ) 
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A4 : 



= -m Hi m Hj (pi| 7 m Pl |p) (pa| 7^ [ps) , 
(Pi\Pr{-4 + m m)Ph \p) (P2\Pr{-4 + m Hi )P L |p 3 ) 



(3.48) 



(Pi|7^b)(p 2 |7^b3) 



(3.49) 



Bl : ( Pl \ ^P L {4+mm)l v PL \p) (p 2 | luP L {-4 + m m ) llx P L \p 3 ) 

= -J (pi\Yl a YPL \P) (jP^laluPL \P3) 

= -V (Pl| 7^L |p> (p 2 | 1,Pl \P 3 ) , (3.50) 
B2 : - (pi\^P L {4 + rn m )P L \p) (p 2 \ Pr(~4 + m m )^P L \p 3 ) 

= -m Hi m Hj 7 M P L \p) (p 2 | 7 M P L |p 3 ) , (3.51) 
B3 : - {p 1 \P R {4 + m m )^P L \p) (p 2 \^P L (-4 + m m )P L \p 3 ) 

= -m Hi m Hj (jp^Pl \p) (p 2 \ I^Pl \ps) , (3.52) 
B4 : ( Pl \ P R (4+ m Hi )P L \p) (p 2 \ P R (~4 + m Hi )P L \p 3 ) 

= - q - ( Pi \YPl \p) (P2\i„Pl \p 3 ) ■ (3.53) 



Thus, all box form factors except P>\ vanish (see Eq. (13. 8j) ). Using the Fierz identity 

(1| YPl |2) (3| 7fl P L |4) = - (3| YP L |2) (1| 7 „P L |4) (3.54) 

and including all the factors, we obtain the generic expressions for the contributions from 
diagrams of types A and B (see Fig. [Hand Appendix [B] for definitions), 



A : B{ 



a 
i 

'An 



E 



- (9%c%!&<£ + <&9%<?li9%) m Hi m Hj D (Ml M 2 , m 2 m , m%) 



(3.55) 



B : B>{ 



a 
> 

'Arc 



E 



49%9ti9'£9£ + j<%c?i<&<% D Q {Ml M 2 , m 2 m , m% 



- {9T2cTA9 3 L2 + <&9%<&.&) m m m Hj D (M 2 , M|, m 2 m , m%) 



(3.56) 



Finally, replacing the vertex coefficients given in Appendix |B] we derive the contributions 
of the heavy gauge bosons and the corresponding would-be-Goldstone bosons: 



P>x{Zh, Z h ) 



a 1 1 v 2 



E^- 



1 + \viVi 1 (.'/,■• </./) - 2 //,■////(■. (.(/,-• '/./) 



1 t; 2 



27rl6s^Af^4/ 



(3.57) 
(3.58) 
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Bi(A H , A H ) 
B^(Zh,A h ) 
with 



a 1 1 1 v 2 
2k 16s^ 25a M w If 



^Xij -3d (y 



a 12 1V 



2tt16 S 4 5M 2 4/2 



^ Xi i \ 3 ^ 



Wa,Vi,yj, 



(3.59) 
(3.60) 



T/ie*T/«M |T/J e |2 
17 /« Vff^l • 



(3.61) 



The SM contribution from the exchange of the light neutrinos is again similar to that from 
Wh, but performing the corresponding replacements. 



Branching ratio 



Collecting everything, the non-vanishing contributions to the vertex and box form factors 
in Eq. (13.1 Op from 7-penguins, Z-penguins and box diagrams in the LHT can be written 



l n 



Q 2 



E V ^ V m 



j^R _ * r M 



Fll 
Flr 



' Mf 



aw 1 v 

a w 1 v 2 

a w 1 - 2s 2 w 1 v 2 
87 2s 2 w M W J 2 

a w 1 v 2 



G w \y i ) + G^(y i ) + -G ( z 1 >(ay i ) 



F w ( 2 / i )+F z ( 2/i ) + -F z (a|/ i ) 
5 



8nM w f 2 
aw I I v 2 



, (3.62) 
(3.63) 
(3.64) 
(3.65) 



in L \ / 

(3.66) 



3 ~ 3 3 ~ 



The branching ratio reads 



B{pL - eee) = Us^MU \A[\ 2 - 2{A\Af + h.c.) + |A?| 2 - f 

^ \ 3 m e 3 

+ i^fl 2 + 1 -{A L X B[* + h.c.) - \(A%B? + h.c.) + i {2\F LL \ 2 + \F LR \ 2 ) 
+ \ {B[F* LL + 2A\Fl L + A\F* LR - AA^F[ L - 2A*F* LR + h.c.) }. 



(3.67) 
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4 Numerical results 



In order to study the bounds on the new parameters imposed by the experimental limits 
on fj, — > e7 and /x — > eee, it is convenient to restrict ourselves to the case of two generations. 
Hence, we are left with four parameters: the LH order parameter /, the masses of the two 
heavy lepton doublets in ( 1A.18I) mm (i = 1, 2), and the angle 9 defining the 2x2 mixing 
matrix between the heavy and the SM charged leptons 

V m = ( ^ V f) = ( C0S9 ). (4.1) 

(In the contributions we consider the e and /i phases, as well as the heavy lepton doublet 
phases, can be safely redefined.) We shall replace m H i, m H2 by 5 and y, however, to 
present our results. The former, which is proportional to the heavy lepton mass difference, 
describes together with 9 the alignment between heavy and SM charged leptons, 

5 = m2 ^- m «\ (4.2) 
m m m H2 

Whereas the latter, which sets the heavy lepton scale, is relevant for discussing decoupling, 

2 

y = Vvm, y i = ^§ L , i = i,2. (4.3) 

Note that both Mw H and mm are proportional to / (see Appendix [A]) . The penguin 
contributions then take the form 

E V m v m HVi) = — [F{vi) ~ nv*)) > (4-4) 

where F stands for a generic function; and the box contributions 
2 sin 26 



E V mKMT F(y t , yj ) = a —{cos 2 9iF( yi , yi ) - F(y 2 , yi )} 

+ sm 2 6[F( yi , y2 )-F( y2 , y2 )}}. (4.5) 



2 



Thus, the LFV amplitudes vanish for vanishing mixing, 9 = 0, or heavy mass splitting, 
5 = 0. 

We plot for illustration in Fig. [5] the form factors for the \x — > and fi — > eee decay 
amplitudes calculated in the previous section as a function of 5 for several 9 and y values and 
/ = 1 TeV. They grow with y and scale like f~ 2 . In contrast with the MSSM case [32,35], 
box contributions to \i — > eee are of the same order than penguins, in particular for y > 1, 
which explains the different behaviour of the decay rates with the sign of 5 for non-maximal 
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Figure 6: Contours of £>(/i — > ej) = 1.2 x 10 11 in the (sin 20, 5) plane for y = 0.25, 1,4 
(left, center, right) and / = 0.5, 1, 2, 3, 4 TeV (from bottom up). 

flavour mixing. The dependence on the new parameters is more clearly seen in Figs. [6] and 
[7J They show the present exclusion contours in the (sin 28, 5) plane implied by the present 
limits on B(fi — > e^y) < 1.2 x 10~ n [38] and £>(/i — > eee) < 10~ 12 [39], respectively, and for 
three values of y, 0.25, 1, 4. The regions above each line of constant / are excluded. As 
it can be observed, mixing angle and mass splitting are correlated, because the alignment 
between the Yukawa couplings of the heavy and the SM charged leptons goes to zero with 
any of them. Present limits on LFV muon decays imply that 9 or 5 < 0.1 for y — 1 and 
/=1 TeV. If no LFV signal is seen by the MEG experiment at PSI, the limits are expected 
to improve by two orders of magnitude [40] and the corresponding exclusion contours would 
be those in Figs. [6] and [7J replacing / by y/lOf. 

As already emphasized, the LFV branching ratios scale like f~ A . However, the y de- 
pendence deserves more discussion. In Fig. [8] we plot the variation of the form factors and 
of the branching ratios with y for maximal mixing, sin 26 = 1, and 5=1. Two comments 
are in order. The non-observation of these LFV processes already sets non-trivial limits 
on the LHT parameters because the central region y ~ 1 is already excluded for natural 
values of the other parameters. More interestingly, B{fi — > eee) goes like y 2 for very large 
y. This is so because y is quadratic in the heavy Yukawa coupling k, which goes to infinity 
with the heavy lepton masses for fixed /. This behaviour is similar to the leading EWPD 
dependence on the top quark mass [41], which scales with m 2 in the region of physical 
interest, allowing a determination of the top mass from a global fit [42]. Just like in the 
top quark case, the dependence is moderate when the particles within multiplets become 
degenerate (the symmetry is recovered). Generic limits from all these figures are tabulated 
in the summary below. 

Let us, finally, comment on the general case with three families. Similarly to the two 
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B(n -> eee) < 10" 




Figure 7: Contours of B(/j, — > eee) = 10 12 in the (sin 29,5) plane for y = 0.25, 1,4 (left, 
center, right) and / = 0.5, 1, 2, 3, 4 TeV (from bottom up). 



20 


xl0~ 6 




15 




Bljl 


10 




/F L r 


5 




A 2 







R 


-5 




A 1 

L 


-10 




\Fll 


-15 






-20 







10 



-2 1Q -1 



10° 

V 



w 1 



w 2 



io- 8 
io- 9 

10 -io 

io- 11 
io- 12 
io- 13 

1Q -14 



sin 26 


1 , iii, , , 

= 5=1 


// = 1 TeV \ 






^e7) - Y 




i . . . .i . . 


> eee) - 



10" 3 10" 2 10" 1 10° 10 1 10 2 10 3 10 4 

y 



Figure 8: Form factors multiplied by (left) and branching ratios (right) as a function 
of y for sin 29 = 5=1 and / = 1 TeV. The latter must be compared with present 
experimental limits on B(fi — > e7) < 1.2 x 10 -11 and B{}i — > eee) < 10 -12 . 
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family case, we have to align the new contributions to the electron and to the muon at the 
10% level. However, now this alignment is not easily related to the usual parameterization 
of the mixing in terms of two mass splittings, three mixing angles and one phase (as in 
the two family case we can safely redefine the e and fi, as well as the heavy lepton doublet 
phases, and then use the same parameterization as for the CKM matrix [42]). In order to 
estimate the fine tuning required by \i —>■ ej, for instance, we rather introduce the ratio 
(see Eq. ( jOD ) 



(4.6) 





3 


2 










i=i 





3 \ / 3 



Era 2 1^)1 Era 2 



.1=1 / \i=l 



which approximately scales like sin 2 295 2 for two families. This is < 10~ 2 for almost all 
the three family parameter space when £>(/i — > ej) is below the present experimental 
limit, and at most ~ 10~ 4 if the limit improves by two orders of magnitude. There is a 
special region in parameter space, however, where the ratio (14. 6 p can be larger even though 
B(fi —>■ 67) is well below the experimental limit. This is around yi ~ 0.3, where the total 
amplitude F(yi) ~ Fw{Vi) + Fz{yi) + \Fz{ayi) in Eq. (I3.32|) is negligible (see right panel 
of Fig. |SJ). But this region is excluded by the present limit on B(fi — ► eee). Thus, the new 
contributions to the electron and to the muon must be aligned at the 10% level, being the 
square of this precision the largest value of the ratio in Eq. (14.61) . Analogously, we can 
define the corresponding ratio using the amplitudes for \i — > eee in Eqs. (I4.4ll4.5p . obtaining 
similar results. The numerical analysis presented here is at some extent complementary 
to the study in Ref. [13], where the correlation between different observables, in particular 
between B([i — *> e^y) and B(fi — > eee), is explicitly shown. 



5 Conclusions 

LH models provide a natural explanation of the little hierarchy between the EW scale 
and the scale where we expect the NP to be, and which is to be explored by the LHC. 
However, these models where the Higgs is a pseudo-Goldstone boson of an approximate 
global symmetry in general suffer some tension in accommodating the many new particles 
required near the TeV scale without upsetting the EWPD constraints. This is ameliorated 
by further extending the model to include a discrete symmetry, the T parity, under which 
all observed particles, including the Higgs boson, are even and hopefully all the new ones 
are odd. All these models, as any universal NP near the TeV scale, must also guarantee 
that the new particles do not mediate too large FCNC processes. We have recalculated 
the new contributions to the LFV processes fi — > ej and [i — > eee in the LHT model, the 
most economical of such proposals. The full Lagrangian has been introduced and all pieces 
of the calculation, in particular the Z-penguin and box diagrams contributing to fi — > eee, 
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B(jjl- 


-> &y) < 1.2 x 10" 11 (sin 


26 = 1) 


f [TeV] 


y = 0.25 


j/=l 


y = 4 


0.5 


\5\ < 0.13 


\5\ < 0.040 


|5| < 0.027 


1.0 


\5\ < 0.52 


|5| < 0.16 


|5| < 0.11 


2.0 


\5\ < 2.2 


|<J| < 0.66 


|<f| < 0.43 


4.0 


\5\ < 14 


\8\ < 3.5 


\5\ < 2.0 



Table 1: Bounds on the splitting 5 from the present experimental limit on £>(/i —>■ ej) for 
several scales / and ratios y, taking sin 26 = 1. 

have been considered in detail. We have found that the former are ultraviolet finite when 
all the Goldstone boson interactions to the order considered are included. Whereas we 
recover previous results for /x — > [12, 13]. 

The present limits on the rates of LFV processes translate into bounds on the LHT 
parameters. Tables [1] and [2] show the bounds imposed in the two family case by /x — > e7 and 
/x — > eee, respectively, on the heavy lepton mass splitting 5 for a maximal mixing angle, 
sin 26 = 1, and several values of the LH scale / and the ratio y related to the common heavy 
lepton mass. The main conclusion is that the new parameters must be tuned to 10% for a 
natural value / ~ 1 TeV. Obviously, raising / quickly reduces the decay rates, which scale 
as f~ A . The results are also sensitive to the parameter y, but the dependence is mild for 
moderate values (see Fig. [SJ, when rates scale roughly like sin 2 26 5 2 . The non-observation 
of LFV effects may be also the result of a conspiracy among the new parameters being 
all slightly above or below their expected natural values, of order one. Analogous fine 
tuning on the alignment of light and heavy leptons is required in the general case with 
three families. 



B{fj. -> eee) < 10" 12 (sin 26 


= 1) 


f [TeV] 


y = 0.25 


y=l 


y = 4 


0.5 


\S\ < 0.16 


\5\ < 0.045 


\5\ < 0.015 


1.0 


\S\ < 0.64 


\8\ < 0.18 


\5\ < 0.061 


2.0 


\5\ < 2.7 


\8\ < 0.72 


\8\ < 0.24 


4.0 


\S\ < 13 


\5\ < 3.3 


\5\ < 0.98 



Table 2: Bounds on the splitting 5 from the present experimental limits on z3(/x — > eee) for 
several scales / and ratios y, taking sin 26 = 1. 

In Table [3] we give both the present and future bounds if the current limits on /x — ► e^y 
and /x — > eee are improved by two orders of magnitude [40]. An asterisk indicates that 
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B(fi -> e 7 ) < 1.2 x 10- 11 (10~ 13 ) 


B(ji -> eee) < 10" 12 (10" 14 ) 


f/TeV > 
sin 26 < 
\6\ < 

y < 


2.5 (8.1) 
0.16 (0.015) 
0.16 (0.015) 
0.16 (*) 


2.3 (7.4) 
0.16 (0.016) 
0.18 (0.018) 



Table 3: Bounds from current (future) experiments on individual LHT parameters, assum- 
ing the others fixed to the natural values / = 1 Te V, sin 26 = 5 = y = 1 . An asterisk 
means that the quoted limit on the branching ratio excludes any y value for the assumed 
values of the other parameters. 

the assumed values are excluded for any possible y. A non-empty region for y is recovered 
increasing / or decreasing sin 26 and/or 8 (see Figs. Eland [7]). Finally, we must note that 
the limits on the corresponding tau decays r — > /r/, ej and r — > fififi, eee are weaker [43], 
typically < 10~ 8 — 10~ 7 . Then, they do not further restrict the order parameter / for 
natural values of the other heavy lepton parameters, but could eventually constrain the 
corresponding mixing angles and heavy lepton masses, which are in principle independent 
of the parameters otherwise involved in the muon to electron processes. However, present 
limits give no significative bound on the parameters related to the third lepton family. 

Note added: 

During the completion of this manuscript several related papers were released. The 
one-loop contributions in the LHT to the tbW vertex in Ref. [44] and to Zll' in Ref. [45] 
have been calculated. In neither case has the order v 2 j f 2 correction to the SM weak boson 
coupling to heavy right-handed fermions been included. More recently a new analysis of 
B decays in the this model has been carried out in Ref. [46], yielding an ultraviolet finite 
result when this correction was taken into account following Ref. [15] and in agreement 
with our findings. 
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A Physical fields 



After the electroweak symmetry breaking (EWSB) the SM gauge boson mass eigenstates, 
which are the T-even, write 



W ± = -^{W 1 TiW 2 ), 
V2 



/ \ / 

Z 



Cw s w 
—s w c w 



\ ( ,\ 

W 3 



\ B I 



(A.l) 



with 



y/2 ' V2 ' 

whereas the T-odd combinations expanding the heavy sector read to order v 2 / f 2 are 



(A.2) 



W^-^iWhTiW 2 ), 



\ Ah J 



( 



1 -XHJj 
V 2 



v 2 \ ( ___„ \ 



(A.3) 



with 



_ W? - Wj _ B 1 -B 2 _ 
w h — 7= ) a n — 7= — > x h — 



5gg' 



V2 



V2 



4(5^2 _ gay 



(A.4) 



Their masses to order v / f are 



M w = ^—\\- 



12f 



, M z = M w /cw, e = gsw = g'cw, v ~ 246 GeV, 



M WH = M ZH =gf[l-^y M AH = g ^(l-^ 2] . 



(A.5) 



The scalar fields must be also rotated into the physical fields [20]: 



7T° 


- 7T° 








- 7T ± 


h 






-> $° 




- $ P 



ir 



12/ 2 
12f 



i2/ 2 y ' 



12/2 



(A.6) 

(A.7) 
(A.8) 
(A.9) 

(A.10) 
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V 



UJ 







$ ± 1 + 



24/ s 



±iu; d 



12/ 



2 : 



5^77 - 4 v / %'(cj° + V2$ p ) - Qgx H uJ } v 2 



24g> 



P 



UJ 



5#(^° + 4 V / 2$ P ) - 4^(5^ + 6g'x H ) v 2 



1 + 



120(7 

.v 2 



P 



, ±i$ ± — 
24/ V P 



(All) 
(A.12) 
(A.13) 

(A.14) 

(A.15) 



For each SM left-handed lepton doublet there is an extra vector- like doublet, 



'it 



( \ 



z = l,2, I 



HR 



Then the left-handed mass eigenstates are 



1L 



I 



2L 



V2 



I 



llL + l 



2L 



HL 



V2 



( \ 



I 



(A.16) 



(A.17) 



where we omit the flavour index. are the SM (T-even) left-handed leptons, whereas 

vhlAhl (vhr,£hr) are T-odd left (right) handed leptons with masses of O(f). The 
SM right-handed fermions are assumed to be singlets under the non-abelian symmetries. 
Heavy leptons receive their masses from the Yukawa term proportional to k f)2.19p . which 
is in general a non-diagonal matrix in flavour space that induces flavour mixing in the T- 
odd sector. The misalignment between the mass matrices of the T-even (SM) and T-odd 
(heavy) sectors is a source of intergeneration mixing (see Section I2T21 . The diagonalisation 
of the k matrix (see Eq. (12.301) ) yield the heavy lepton masses 



rrto 



m m [l- — 



{k.li 



B Feynman rules 

We present below just the Feynman rules which are necessary for the calculation of the 
LFV processes discussed in this work. They are given in terms of generic couplings for the 
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following general vertices involving scalars (S), fermions (F) and/or gauge bosons (V): 



[V,FF] 


= ie^(g L P L + g R P R ), 


(B.l) 


[SFF] 


= \e{c L P L + c R P R ), 


(B.2) 




= ieKg^, 


(B.3) 


[V M S( Pl )S( P2 )] 


= ieG(pi - p 2 f, 


(B.4) 


[V M (p 1 )V,(p 2 )V p (p 3 )] 


= ieJ [g^(p 2 - Pi y + g^(p 3 - p 2 y + g^( Pl - p 3 )1 


(B.5) 



where all momenta are assumed incoming. The conjugate vertices are obtained replacing: 
9L,R^g* L , R , c LjR ^c RL , K <-> K*, G^G*, J^J*. (B.6) 

B.l SM with massive neutrinos 



For comparison we first give the rules for the SM with light massive neutrinos. The sign 
conventions for the covariant derivatives are those in Ref. [22]. 



VFF 


ifP 


zfp 






9l 
9r 


-QfSij 
-QfSij 


Z L 5ij 
ZftSij 


/X PMNS 

V 2s w 



1 v ji 

fK PMNS 

V 2sw 




where Z[ R = (vf ± df)/2swcw with t>/ = T( l — 2<5/s^ and a/ = T/ J 



SFF 



! TTl^s -irji* 

"7^^ PMNS 

I ^P_Trji* 



svv 






K 







VSS 






G 




r 2 - s 2 
IswCw 
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vvv 


7 w + w 


zw+w~ 


J 


-1 


c-w/ sw 



The fields ± are the would-be Goldstone bosons eaten by the gauge bosons fields W 
after the EWSB. 



B.2 LHT model 



The sign conventions are chosen to be compatible with those employed for the SM (which 
coincide with those in [26] up to a sign in the definition of the abelian gauge couplings 
in the covariant derivative in Eq. (12.141) ). In particular, these Feynman rules include the 
0(v 2 /f 2 ) contribution to the Zv % H v J H vertex missed in the literature. 



VFF 



9l 
9r 



VFF 



9l 
9r 



ifhfu 



-Q/Sij 
-QfSij 



Zv^ H 



— Sij 

2s\yC\y 

-J-fi-Hl 

2s w c w \ Ap 



'10 



zt H e H 



i 



1 

2s\vCyy 



(-l + 2s 2 w )Si 
-l + 2s 2 w )5i 



A 



lOcw 



2s w P 




2s w ' lOcwP) m 




1 



V2s w 




SFF 


T]t H & 


0J°t H P 


cl 
cr 


i m t H 


' v 2 (5 s w } 

i m e i yij 
10c w M Ah m 




yij 
V H£ 


i m e H 


" , v 2 ( 1 c w \ 

i rn.fi ij 
2s w M Zu m 




yij 


10c w M Ah 


2s w M Zh 



SFF 




cl 
cr 


i m v i 


V2sw Mw H 
i rn.fi ij 

\/2sw M\v H H 
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svv 






K 


±iM WH 





vss 






G 




± cw 
s w 


K 1 8 pel) 



VVV 




ZW H W H 


J 


-1 


c-w/ sw 



The fields u ± , u° and r\ are the Goldstone bosons of the [SU(2) xU(l)]i x [SU(2) xU(l)] 2 
breaking into its diagonal subgroup. They are eaten by the heavy gauge bosons W H , Z H 
and Afj, respectively. (Actually these Goldstone bosons mix with an additional physical 
Higgs triplet $ at order v 2 /f 2 and it is this linear combination of fields that is eaten.) 
In principle, also the scalar triplet $ contributes to the processes considered here. The 
corresponding diagrams can be obtained replacing by $ ± and Zh, Ah by <3>° and $p. 
The Feynman rules for the vertices containing $, neglecting the masses of the SM fermions, 
involve couplings of 0{y 2 / f 2 ). As each diagram contains at least two such vertices, if any, 
they are suppressed by a factor of 0(v 4 / f 4 ) [26]. 

C Loop integrals 



Pl P2 




Figure 9: Generic one-loop diagram with iV legs. 
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Consider the generic one-loop diagram with N legs in Fig. [9J where 



N-l 



h=Pi, h=Pi+P2, ■■■ fcjy-i = /,Pi- (C.l) 



i=l 



This diagram involves integrals of the type 

' rpN _ A-D f %l " " " Qnp /p r,\ 



167r 2 mi -mp ^ y (27r) D (g 2 -mg)[(g + Jfei) 2 -m?]---[(g + A;iv-i) 2 -Tn 



iV- 



These integrals are symmetric under permutation of the Lorentz indices. The integration 
is performed in dimensional regularization. The mass scale \i keeps track of the correct 
dimension of the integral in D = 4 — e spacetime dimensions. P < N is the number of 
g's in the numerator and determines the tensor structure of the integral (scalar for P = 0, 
vector for P — 1, etc.) The notation is A for T , B for T 2 , etc. and the scalar integrals are 
Aq, Bq, etc. The tensor integrals can be decomposed into a linear combination of Lorentz 
covariant tensors constructed from g^ v and a linearly independent set of the momenta [47] . 
The choice of the basis is not unique. Here we choose g^ v and the momenta k{, which are 
sums of the external momenta Pi [22]. In this basis, the tensor-coefficient functions are 
totally symmetric in their indices. For this work, we need the following decompositions: 

= k^B, , (C.3) 

c, = k^a + k^Ct, (c.4) 
2 



Cy U — dftvCoo + ^ kiykj U Cij , (C-5) 

3 

D„ = ^h^Di, (C.6) 



Dyu — g^uP J oo + kiykj U Dij . (C.7) 

These functions have been calculated for the argument configuration required by the 
processes under study, obtaining the following results. 



C.l Two-point functions 

Consider now the diagram with two legs in Fig. [TUJ 

{1, <," 1 



i f d q 

,^{4 B^} (args) = / 

I67H y (27r)^ (g J — TOq) [(g + p) J — m\ 



21 ' 



(C.8) 



35 




Figure 10: Generic one-loop diagram with two legs. 

where k\ = p. The corresponding tensor coefficients are functions of the invariant quantities 
(args) = (p 2 , ml, m\). The functions B = B(0; M 2 , Mf) and B = 5(0; M 2 , Mf) read 



B = B C 



A 6 + l- 



M'( In ^ - M| In 



2^ Ml 



fJL' 



Mf - M.j 



(C.9) 



2 



M 2 M 2 
4M 2 Mf - 3Mf - M 2 4 + 2M 4 In — + 2M 2 (M 2 - 2M 2 ) In ■ 2 



4(M 2 - M 2 ) 2 



with A f 



: -B -B 1 , (CIO) 
7 + ln47r. These functions are ultraviolet divergent in D = 4 dimensions. 



C.2 Three-point functions 




P2 - Pi 



"P2 



Figure 11: Generic one- loop diagram with three legs. 
Consider now the diagram with three legs in Fig. [TO 



16n 



2 {C , C, C^} (args) 



4--D 



d D g 



{1, g", g>V} 



(2tt) d (g 2 - mg) [(g + pi) 2 - mf] [(g + p 2 ) 2 - m|] 



21 ' 



(C.ll) 
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where we have chosen the extenal momenta so that k\ — pi, k 2 — Vi- The corresponding 
tensor coefficients depend on the invariant quantities (args) = (p 2 , Q 2 -,p\\ m 2 ), m 2 , m?), with 
Q 2 = (p 2 - Pl ) 2 . The functions C = C(0, Q 2 , 0; Ml, Mf, Mf) with x = M 2 /M 2 read 



C 



1 

M 2 



x + In a; 



(1 -xf 

Q 2 -2 - 3a; + Qx 2 - x 3 - 6xlnx 



Cn = C*22 = 2 C12 

Coo 



M 2 12x(l - x) 4 

1 -3 + Ax - x 2 - 21nx 



+ 0(Q 2 ), 



M 2 4(1 - x) 3 
1 11 - 18x + 9x 2 -2x 3 + 61nx 



M 2 



18(1 -x) 4 
Q 2 ll-18x + 9x 2 -2x 3 + 61nx 



Ml 



72(1 -x) 4 



+ 0(Q 4 ), (C.12) 
(C.13) 

+ 0(Q 2 ), (C.14) 
+ G(Q 4 ). (C.15) 



Or else, defining C = C(0, Q 2 , 0; M|, M 2 , M 2 ), 

-1 + a; — x In x 



Co 



1 

Ml 



(1-x) 2 

Q 2 -1 + 6x - 3x 2 - 2x 3 + Qx 2 lnx 



Cl = C*2 
Cn = C*22 = 2 C12 

Coo 



M 2 12(1 - x) 4 

1 1 - 4x + 3x 2 - 2x 2 In x 
Mf 4(1 -x) 3 ' 
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+ 0(Q 4 ), (C.16) 



Ml' 



2 1 



18(1 -x) 4 
Q 2 -2 + 9x - 18x 2 + llx 3 - 6x 3 lnx 



Ml 



72(1 -x) 4 



(C.17) 
(C.18) 
+ 0(Q 4 ).(C.19) 



Note that C 00 and C o are ultraviolet divergent in D = 4 dimensions. 

In the limit Q 2 = the following useful relations among two- and three-point functions 



hold: 



B 1 + 2C 



oo 



X 



± + ±B 1 + C 00 -^M-fC 



0, 
0, 



1 



+ B l + 6C 00 - xMlCo = A e -ln 



T 2 



(C.20) 
(C.21) 

(C.22) 
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C.3 Four-point functions 

The ones we need are all ultraviolet finite: 

{A,, £>", £>""} (args) = 
d 4 g 



16tt 2 



{1, q», q»q»} 



(2tt) 4 (g 2 - m 2 ) {{q + k x f - m\} [(q + k 2 ) 2 - m 2 ] [{q + k 3 ) 2 - 



m 



21 > 



(C.23) 



3 

with kj = ^Pi and (args) = (pj, pi, pi, pi, (pi +P2) 2 , (P2 + P3) 2 ; ml, m 2 , m\, m 2 ). In the 

i=l 

limit of zero external momenta, only the following integrals are relevant: 



16^° " 
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£>nn = - 



d 4 g 
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(C.24) 
(C.25) 



In terms of the mass ratios x = m\jm\, y = m\j 'ml, z = m\jm\ the integrals above can 
be written as: 



d (x, y, z) = m^Do = 



xlnx 



ylny 
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For two equal masses (mo = m 3 ) we get 
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